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We present the first measurement of the color representation of the hadronically decaying W 
boson in ti events, from 5.3 fb^^ of integrated luminosity collected with the DO experiment. A novel 
calorimeter-based vectorial variable, "jet pull," is used, sensitive to the color-flow structure of the 
final state. We find that the fraction of uncolored W bosons is 0.56 ± 0.42(stat-|-syst), in agreement 
with the standard model. 

PACS numbers: 12.38.Qk, 12.38.Aw, 14.65. Ha 



Color charge is conserved in quantum chromo dynamics 
(QCD), the theory that describes strong interactions 
At leading order in the strong coupling constant as, color 
can be traced from initial partons to final-state partons 
in high-energy hadron collisions. Two final-state partons 
on the same color-flow line are "color-connected" and 
attracted by the strong force. As these colored states 
hadronize, the potential energy of the strong force be- 
tween them is released in the form of hadrons. Thus, 
knowledge of the color-connections between jets can serve 
as a powerful tool for separating processes that otherwise 
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appear similar. For example, in the decay of a Higgs {H) 
boson to a pair of bottom (6) quarks, the two h quarks are 
color connected to each other, since the H is uncolored 
(color singlet), whereas in g — )■ 66 background events, 
they are color-connected to beam remnants because the 
gluon carries a color and an anti-color (color-octet). We 
follow a recent suggestion Q for reconstructing these 
color connections experimentally, using observables that 
can be modeled reliably by available leading-log parton- 
shower simulations. The technique involves measuring a 
vectorial quantity called "jet pull," which represents the 
eccentricity of the jet in the rj-(j) plane Q and the di- 
rection of the major axis of the ellipse formed from the 
jet energy pattern. Jets tend to have their pull pointing 
towards their color-connected partner. For instance, in 
H ^ bb events, the pulls of the two 6-jets tend to point 
towards each other, whereas in g ^ bb events, they point 
in opposite directions along the collision axis. 

Verification of color flow simulation and jet pull recon- 
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FIG. 1: (Color online) Diagram showing two jets in the rj-ff) 
plane, and the reconstruction of the jet pull vectors (t), jet 
pull angles (6)^""), and relative jet pull angles (6)^""). 



struction for both color-singlet and color-octet configu- 
rations is interesting in its own right jj] and is needed 
before jet pull can be used in, e.g., H ^ hh searches. 
Color-octet patterns can be studied in many processes, 
such as W / Z boson production in association with jets. 
A pure sample of color-singlet hadronic decays is diffi- 
cult to obtain at a hadron collider, but ti events with 
an £-|-jets final state are good candidates since they have 
a characteristic signature and contain two jets from the 
decay of a boson, which is a color singlet. Each of 
the two 6-jets coming from the top quark decays is color- 
connected to one of the beam remnants in a color-octet 
pattern. 

In this Letter, we use data collected with the DO detec- 
tor at the Fermilab Tevatron pp collider, correspond- 
ing to 5.3 fb^^ of integrated luminosity, to present the 
first experimental results on the study of jet pull, using ti 
events decaying to £-|-jets [ti WhWb — >■ ivbjjb, where 
£ — e,fi). The object identification, event selection, and 
simulated Monte Carlo (MC) events are the same as those 
used in the tt cross section analysis [6\ , except that looser 
6-tagging criteria are used to increase the statistics of 
double 6-tagged events. We obtain a ~ 90% pure tt sam- 
ple by requiring an isolated lepton with px > 20 GeV, 
missing transverse energy 1^t> 20 GeV (> 25 GeV for 
the /i-|-jets channel), and at least four jets, reconstructed 
with a midpoint cone algorithm 8] of radius 0.5, with 
Pt > 20 GeV. At least one jet must have pr > 40 GeV, 
and at least two jets must be identified as 6-jets. Table |T] 
shows the event yields for these selection criteria. 

To extract the fraction of color-singlet hadronic W bo- 
son decays, the data are compared to both standard 
model tt MC (with a color-singlet W boson) and an 
alternative model of tt with a hypothetical color-octet 
'W boson decaying hadronically with identical proper- 
ties except for its color representation. The latter is sim- 
ulated using the MadGraph (MG) 01 event generator 
interfaced to PYTHIA for showering and hadroniza- 
tion. Simulated events are processed with a GeantS- 



TABLE I: Yields of events passing selections with exactly 4 or 
> 5 jets. At least two 6-tagged jets are required in the anal- 
ysis, but the numbers of events with zero or one 6-tagged jet 
are also given. The number of ti events is calculated using the 
cross section determined with this data sample, = 8.50 pb. 
Uncertainties include statistical and systematic contributions. 
The total uncertainties are smaller than the sum of individual 
uncertainties due to negative correlations between samples. 



channel 


sample 


6-tags 


1 6-tag 


> 2 6-tags 


^-l-4jets 


W^-l-jets 


576 ± 75 


229 ± 32 


49 ± 8 




Muhijet 


115 ± 16 


46 ± 7 


7 ± 2 




Z-l-jets 


42 ± 6 


16 ± 3 


4 ± 1 




Other 


31 ± 4 


19 ± 2 


9 ± 1 




ti 


160 ± 22 


417 ± 38 


519 ± 51 




Total 


923 ± 62 


727 ± 24 


589 ± 48 




Observed 


923 


743 


572 


^+>5jets 


W^-l-jets 


60 ± 22 


26 ± 11 


7 ± 3 




Multijet 


17 ± 3 


12 ± 2 


3 ± 1 




Z-l-jets 


4 ± 1 


2 ± 1 


1 ± 1 




Other 


3 ± 1 


3 ± 1 


2 ± 1 




tt 


34 ± 6 


90 ± 13 


132 ± 17 




Total 


118 ± 19 


132 ± 7 


145 ± 15 




Observed 


112 


127 


156 



based \lv\ detector simulation, overlaid with random data 
to account for backgrounds, and reconstructed as data. 

DO uses three liquid- argon/uranium calorimeters to 
measure the energies of particles: a central section (CC) 
covering \ri\ up to « 1.1 and two end calorimeters (EC) 
that extend coverage to \r]\ « 4.2 housed in separate 
cryostats In addition, scintillators between the CC 
and EC cryostats provide sampling of developing show- 
ers for 1.1 < |?7| < 1.4. There are approximately ten 
layers in the radial direction (depending on ry), generally 
composed of cells spanning 0.1 x 0.1 in t] x (j>. The en- 
ergy resolution is about 15%/V^ ® 0.3% (in GeV) for 
electrons and 50%/VE (B 5% for hadrons. Pileup energy 
from overlapping pp interactions result in about 0.5% of 
cells having energy above the noise-limited energy thresh- 
old (« 50 — 500 MeV, depending on layer and ry). This 
energy is roughly exponentially distributed, with a mean 
of « 350 MeV. 

The pull is determined for each jet of a pair of re- 
constructed jets, using the measured energies of the 
calorimeter cells (see Fig. [T]). Each cell within AR — 
y/ {Acj))'^ + (ATy)^ < 0.7 of the i?T-weighted center of one 
of the jets of the pair {rf^^ , is assigned to the jet 

nearer in AR. The contribution of each selected cell 
to the jet pull is iccii = E^^^fcci\\rccih where fccW = 
ivT^^ - Vd^, 0™" - 'A'"'), and Eij?^^ is the ceh's transverse 
energy with respect to the nominal center of the detec- 
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FIG. 
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2; (Color online) The average of the two jet 6'^^^ dis- 
tributions for jets in pairing (a) w and (b) b, in events with 
exactly four jets, at least two 6-tags, and the Mw requirement 
on the ui-pair jets. The x^/ndf compares the data to the total 
MC distribution. 
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FIG. 3: (Color online) (a) Leading-pT and (b) second-leading- 
Pt jet distributions for «;-pair jets, in events with two 

jets and no 6-tagged jets. The x^/ndf compares the data to 
the total MC distribution. 



tor. The jet pull ist — cells i ^/-^t*- "^^^ polar angle of 
the jet pull, 0P"", is defined to be zero when pointing in 
the positive 77 direction along the beamline. A small cor- 
rection to the jet pull is made to account for the energy 
response and noise in the calorimeters as a function of rjd, 
particularly in regions between the central and forward 
cryostats. The angle of the jet pull direction relative to 
the line defined by the centers of the jet pair (6'^'^"") is 
also of interest, as we expect color-connected jets to have 
pulls pointing towards each other. The 9^^^^ quantity is 
calculated for each jet in the pair of highest-px 6-tagged 
jets (b pair) and the pair with highest px which are not 
amongst the two highest px 6-tagged jets {w pair). 

To select events with a higher purity of properly identi- 
fied jet pairs from hadronic W boson decays, we split the 
sample into events where the invariant mass of the w-pair 
jets is consistent with the W boson mass, \mjj — Mw\ < 
30 GeV, and events where it is not. For the former, these 
two jets are found to match the partons from the W bo- 
son decay within AR < 0.5 in 66% of tt MC events with 
four jets and 46% of events with 5 or more jets. In the 
latter case, additional gluon radiation in the initial or fi- 
nal state leads to possible additional color configurations, 
diluting the measurement. 

Since the w-pair jets in ti events are often from the 
W boson decay, we expect them to be color-connected, 
thus the jet pulls should generally point towards each 
other. We expect 6-pair jets to have one of the &-jets 
color-connected to the proton beam and the other to the 



anti-proton beam, thus the jet pulls should be generally 
pointing away from each other. This tendency is seen in 
data as shown in Fig. O with smaller 9f^^^ in the w pair 
than in the b pair. However, the jets in w and b pairs have 
different kinematics, separation in the detector, and fla- 
vor. A direct interpretation of the effects from color-flow 
is therefore not possible from this comparison. Further- 
more, there are detector and reconstruction effects on jet 
pulls from overlapping jet pull cones, calorimeter noise 
and pileup, and calorimeter response inhomogeneity. For 
instance, there would be fewer cone overlaps if the jet pull 
was defined using only calorimeter cells within AR < 0.5, 
producing on average smaller values for 0^^^^. With this 
alternative definition the shape in Fig. [2][a) would peak 
more towards zero and that in Fig. [2][b) would be flatter. 
These effects are found to be well-modeled by the simu- 
lation, and the jet pull definition based on the AR < 0.7 
cone gives a slightly improved singiet-octet separation. 
The relative jet pulls 9f^^^ in data are also found to be 
well-modeled by simulation for other jet pairings, such as 
a random w-pair jet and a random 6-pair jet. In control 
samples consisting of events with a leptonic W boson de- 
cay, and two, three, or four jets, none identified as 6-jets, 
various jet pairings also have jet pulls that agree with 
simulations. Figure |3] shows the 9f^^^ distributions for 
jets in a control sample with a leptonic W boson decay 
and two not-fe-tagged jets. 

To quantify the method's sensitivity to the color- 
fiow structure (color-singlet versus color-octet) for the 
hadronic W boson decay, we fit the data to two hypothe- 
ses: (i) standard model tt with a color-singiet hadroni- 
cally decaying W boson (singlet MC) and (ii) tt with a 
hypothetical color-octet 'W boson (octet MC). We de- 
termine the fraction of events coming from color-singlet 
W boson decay (/singict) using the fitting procedure from 
the DO combined tt cross section analysis We simul- 
taneously measure the tt cross section to avoid any possi- 
ble infiuence of the tt signal normalization on the /singlet 
measurement. The discriminating variable used for the 
fit is derived from the 6'J|,"" angles of the w-pair jets and 
depends on the AR between the two jets and their 77^. For 
events failing the W mass requirement, wc do not split 
the regions further; for other events we split the data 
sample according to the r]d of the jets and AR between 
the jets. For events where the two jets are highly sepa- 
rated {AR > 2), we use the 0^^,"" of the leading-pT jet. 
Little discrimination is possible for these events, since the 
additional color radiation is distributed over a large area 
of the calorimeter. When the two jets are close {AR < 2) 
and \rid\ < 1.0 for both jets, we use the minimum 9^^^ of 
the two jets. This is the most sensitive region, and the jet 
pull is accurately reconstructed in the central calorime- 
ter due to less pileup energy and uniformity of response. 
Otherwise, if \?]d\ of the leading-pT jet is < 1.0 (> 1.0), 
the 9^^^^ of the leading-py (second- leading pt) jet is used. 
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TABLE II: The one standard deviation (a) variation of /singict 
from main systematic uncertainties. The total systematic un- 
certainty includes all uncertainties, summed in quadrature. 



Source 


+ 1(7 


-la 


Singlet/octet MC shapes 


0.188 


-0.188 


Jet pull reconstruction 


0.100 


-0.093 


Jet energy resolution 


0.033 


-0.013 


Vertex confirmation 


0.028 


-0.029 


PYTHIA tunes 


0.023 


-0.025 


Jet energy scale 


0.024 


-0.009 


Jet reconstruction and identification 0.017 


-0.017 


tt modeling 


0.014 


-0.033 


Event statistics for matrix method 


0.009 


-0.010 


Other Monte Carlo statistics 


0.009 


-0.007 


Multijet background 


0.006 


-0.007 


Total systematic 


0.222 


-0.218 
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FIG. 4: (Color online) The discriminating color-flow variable, 
the minimum Sf^,"" for the w-pair jets, for events passing the 
Mw requirement, with AR < 2, and rja < 1.0 for both jets. 
The tt MC shape is obtained using the measured value of 

/Singlct. 



Table |TI] lists the contribution of each non-negligible 
source of systematic uncertainty on /smgiot- For all but 
the theoretical cross sections, MC statistics, and normal- 
ization of the heavy flavor jets background uncertain- 
ties, we apply the systematic uncertainties just to the it 
signal sample and ignore the effect on background, as the 
purity of the ti sample is high. To estimate the possible 
systematic shift of the 9f^^^ distribution due to the differ- 
ent energy scale and noise of the calorimeter cells between 
data and MC as a function of rjd, we apply ±50% of the 
jet pull 77 correction and take the resulting difference in 
shape as the systematic uncertainty for jet pull recon- 
struction. This covers the differences in the average 6*^"" 
when comparing data and MC control samples. We also 
study systematic uncertainties as in the main ones 
being from the jet energy scale, jet energy resolution, 
6-tagging efficiency, and lepton misidentification. Addi- 
tional systematic uncertainties on df^i^ are assessed to 
account for possible differences between MC and data re- 
lated to the modeling of underlying event, hadronization, 
and jet showering. To estimate the variation due to these 
possible mis-modelings, we compare 0^^^^ distributions in 
events simulated with pythia to those with ALPGEN [i3| 
or MC@NLO [ill, and showering with herwig 15 |. We 



also do the comparisons for various pythia parameters 
for underlying event and color- reconnection 16|, such as 
tunes APro and NOCR 17|. When deriving /singiot from 
the fit, we use the maximal variation obtained with the 
different 9f^^^ distributions as an estimate of the system- 
atic uncertainty. 

Since the results are statistically limited and the anal- 
ysis does not as yet provide sufficient sensitivity for 
a definitive observation of color-flow, we set limits on 
/singiot using the likelihood ratio ordering scheme of Feld- 



man and Cousins [18| . We follow the same approach used 
for the simultaneous extraction of the ratio of branch- 
ing fractions and the tt cross section [l^ and generate 
ensembles of pseudo-experiments for different values of 
/singiot between and 1, with the tt cross section fixed to 
the measured value. We then vary the systematic uncer- 
tainties using Gaussian distributions and perform the fit 
as for the measurement on data. Statistical uncertainties 
are incorporated by smearing the measured value for each 
pseudo-experiment with the uncertainty determined in 
data. We use the nuisance parameters method where the 
expectation is fit to the data, for a variation of the initial 
prediction within the systematic uncertainties, allowing 
also the central result to change @ . Other methods give 
compatible results. 

We measure /singiet = 0.56 ± 0.42 [±G.36(stat) ± 
0.22(syst)] and att = 8.50to;76 pb, consistent with our 
dedicated cross section measurement . Figure |4] shows 
the distribution for one of the regions of the discriminat- 
ing color-fiow variable, using the measured tt cross sec- 
tion and measured /singiet. The expected 99% C.L. and 
95% C.L. limits are /singiot > 0.011 and /singiot > 0.277 
respectively, corresponding to an expected sensitivity to 
exclude /singict = of about three standard deviations, 
based on pseudo-experiments. The 68% C.L. allowed re- 
gion from data is 0.179 < /singict < 0.879. Figure[S]shows 
the expected 68%, 95%, and 99% C.L. bands for /singict- 

In summary, we have presented the first study of color 
flow in tt events, with the method of jet pull, using 
5.3 fb^^ of DO integrated luminosity. The standard 
model MC predictions are found to be in good agree- 
ment with data, for both the jets from the hadroni- 
cally decaying W boson, which should be in a color- 
singlet configuration, and the 6-tagged jets from the 
top quark decays, which should be in a color-octet con- 
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FIG. 5: (Color online) Expected C.L. bands for /singiet- The 
measured value is shown on the horizontal axis, and the input 
value on the vertical axis. The wide-dashed line shows the 
expected value and the black-white fine-dashed line indicates 
the measured value of /singict. 



figuration. To quantify our ability to separate singlet 
from octet color-flow, we measured the color represen- 
tation of the hadronically decaying W boson and found 
/singict — 0.56±0.42(stat-|-syst), while the expected 95% 
C.L. limit was /singiet > 0.277. The ability to use color 
flow information experimentally will benefit a wide range 
of measurements and searches for new physics. 
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